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Abstract 

An array of planar Penning traps, holding single electrons, can realize an artificial molecule 
suitable for NMR-like quantum information processing. The effective spin-spin coupling is accom- 
plished by applying a magnetic field gradient, combined to the Coulomb interaction acting between 
the charged particles. The system lends itself to scalability, since the same substrate can easily 
accomodate an arbitrary number of traps. Moreover, the coupling strength is tunable and under 
experimental control. Our theoretical predictions take into account a realistic setting, within the 
reach of current technology. 



I. INTRODUCTION 



In this proposal we bring together the best of two avenues to quantum information pro- 
cessing: nuclear magnetic resonance (NMR) Q| and ion trapping j^. Both approaches have 
provided the first experimental demonstrations of fundamental quantum logic gates and 
quantum algorithms, though still limited to few qubits. Actually, it is relatively simple 
to build a small NMR quantum computer, useful for a proof-of-principle test of quantum 
algorithms, but unable to perform any real computation, involving hundreds of qubits. In- 
deed, the scalability seems to be a rather tough question for NMR quantum computing. 
Another disadvantage of NMR quantum computation is that NMR experiments deal with a 
large number of molecules, building up an ensemble of indistinguishable quantum computers. 
This fact brings in relevant theoretical and practical implications: from the system initialisa- 
tion to the debate on the same quantum character of the computation carried out with such 
a device. However, NMR techniques, based on radio frequency (rf) and microwave (mw) 
pulses, enable to carefully prepare, manipulate, and detect the qubit states with relative 
ease. Spatial separation of the qubits is not required, since different qubits are distinguished 
using different resonance frequencies. 

On the other hand, in ion trap quantum computation, qubits are stored in isolated 
quantum systems, arranged to form strings of trapped particles, spatially separated and 
singly addressable with optical radiation Q. Such systems can be prepared in their 
motional ground state via sophisticated laser cooling techniques. Coherent manipulation 
of the qubits requires strongly focused pulses of controlled intensity, phase and duration at 
optical frequencies. The required setup is rather involved and the experimental realization 
challenging. These technical problems have motivated other proposals by Wunderlich et 
al. to implement a quantum computer, based on trapped ions in a linear Paul trap, but 
using long- wavelength radiation, in the radio frequency or microwave range \^\. To this 
end, internal and external degrees of freedom of the trapped ion are coupled by means of a 
magnetic field gradient. Moreover, in a linear Paul trap, the collective vibrational modes of 
a chain of two-level ions extend this coupling to different ions. The system is formally 
analogous to a collection of spins, interacting through the so-called J-coupling, typical of 
nuclear spins in molecules However, the proposal by Wunderlich et al. presents 

some drawbacks, especially for the scalability, due to the fact that all the ions are stored 
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in the same linear trap. A more flexible design would be based on a string of individually 
tailored microtraps, each of them trapping a single ion |^ . These ideas are versatile and can 
be adapted to other scenarios. 

Our system consists of an array of Penning traps, each of them confining a single electron. 
A Penning trap makes use of static electric and magnetic fields to trap charged particles, 
like ions or electrons. In particular, the magnetic field provides the radial confinement, 
which, in a Paul trap, is achieved by means of an oscillating (radio frequency) electric field. 
The resulting dynamics of a particle in a Paul trap is harmonic in all directions, whereas 
in a Penning trap the combination of the electric and magnetic fields gives rise to a more 
complicated orbital motion, which is the superposition of the magnetron and cyclotron 
oscillators. These differences should be taken into account, when describing the effects of 
an additional magnetic gradient on the trapped particle dynamics. 

We choose to trap electrons instead of ions, because of their smaller mass, which results 
in higher trapping frequency for the quantized external degrees of freedom. Indeed, the 
typical resonance frequencies of the resulting electron motion lie in the radiofrequency and 
microwave domain, making it possible to employ the same technological resources and meth- 
ods developed for NMR experiments. Moreover, differently from a Paul trap, a Penning trap 
does not rely on rf fields, a benefit in terms of stability of the trapping potential. In turns, 
this translates into less decoherence affecting the trapped particles. Therefore, electrons 
trapped in vacuum seem to be a promising candidate for quantum information processing 
[10, llll , taking advantage of the techniques and strategies devised both for NMR and ion 
trapping quantum computing. 



What we have in mind is a new concept of planar Penning traps characterised by 
an open geometry. The trap electrodes are deposited on a ceramic substrate by means of 
well-established thin film or thick film technology, which allows for a variety of different con- 
figurations as well as dimensions. The trapping mechanism relies, as in conventional Penning 
traps, on the application of a magnetic field together with an electrostatic quadrupole po- 
tential. A single trapped electron is confined in vacuum at an adjustable distance from the 
trap surface. The same substrate can accomodate several traps in order to form a regular 
one or two dimensional array of trapped particles. Qubits are encoded in the natural two- 
level system provided by the electron spin in the external strong magnetic field, similarly 
to what happens for NMR spin one-half nuclei. The two possible spin orientations | t) and 
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I I) represent, respectively, the logic states |1) and |0). Here, however, the two spin levels 
are greatly separated in energy and thermal excitation is completely negligible, especially 
at the trap cryogenic temperatures. Hence, the system, after initialisation, remains in its 
ground state, corresponding to the spin-down state. We have already observed that in order 
to make the spin qubits distinguishable, one has to introduce a small magnetic gradient Q| . 
Actually, with a judicious choice of the magnetic gradient, that in our configuration depends 
on all the spatial coordinates, we can also build up a molecule-like system. In the present 
geometry, the magnetic field gradient is applied across the substrate in order to differentiate 
among the spin resonance frequencies at each trapping site. Therefore, the spin qubits are 
distinguishable and can be selectively frequency addressed via microwave pulses. The same 
magnetic field gradient, mediated by the long range Coulomb interaction, enables the effec- 
tive spin-spin coupling between different electrons. The resulting system may be regarded 
as an artificial molecule suitable for NMR quantum computation. Actually, we can even 
envisage applications to simulate other quantum systems, like the Ising model. Our system 
offers obvious advantages in terms of scalability. In addition, the spin-spin coupling depends 
on external parameters, like the strength of the magnetic field gradient, the trapping fre- 
quencies, and the trap separation, that can be adjusted to obtain the optimal performance 
of the quantum processor. We point out that this coupling is proportional to 1/cP, with d 
being the inter-trap distance. Therefore, the spin-spin coupling strength is relevant only for 
nearest-neighbor electrons, while rapidly decays along the electron chain. 

Finally, for what concerns the initialisation of the system, or, in other words, how to 
reset the quantum register to the state |0), several strategies could be followed. A possible 
one is to apply a transverse oscillating magnetic field resonant with the frequency difference 
between the cyclotron motion and the spin precession. When the magnetic gradient is off 
the transverse field will flip the spin up state of each particle transfering the energy to 
the corresponding cyclotron motion, which will release its energy to the environment via 
syncrotron radiation. 

The manuscript is organized as follows. In Sec. |H] we describe the trap design and 
how to create a scalable device, putting several planar traps on the same substrate. The 
theoretical framework is developed in Sec. IIIH where we discuss the role of the applied 
magnetic field gradient to achieve the individual frequency addressability of the qubits (Sec. 
nil A|l and derive an analytical expression for the effective spin-spin coupling (Sec. IIIIBj) . As 
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an example we illustrate how to implement a set of fundamental quantum logic gates with 
the resulting NMR-like molecule, made out of trapped electrons (Sec. IIV|) . The concluding 
remarks and the future perspectives for the system are summarized in Sec. |V1 



II. THE PLANAR PENNING TRAP 

We have tried to design a trap which fulfills the needs of quantum computation and at the 
same time ensures the same amount of control and precision already achieved with conven- 
tional traps. To that end we consider a planar Penning trap, which is shown schematically 
in Fig. ^ The planar trap is a Penning trap in the sense that both electric and magnetic 
static fields are used to confine three-dimensionally a charged object. A generic Penning 
trap generates a quadrupolar static electric field together with a homogeneous and static 
magnetic field aligned in the direction of the symmetry axis, which we'll assume to be z. The 
magnetic field confines a c harg ed particle radially, whereas the electric field provides the con- 
finement along the z axis Since the ideal quadrupole potential depends on the square 
of the coordinates, the motion in an ideal trap can be decomposed in three independent 
harmonic oscillators which are commonly denominated by cyclotron, magnetron and axial, 
each of them having different characteristic frequencies. When the trap is not ideal, due to 
imperfectly homogeneous or quadrupolar fields, respectively, the three eigenmotions couple 
and can no longer be described as independent. However, the effects of such a coupling can 



be used for a number of applications, see e.g. (l4l. llSl Il6|. 
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FIG. 1: Simplest configuration of a planar trap with two electrodes (black) on an isolator substrate 
(white) . 
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A. General trap properties 

A planar trap is a novel concept of trap which allows for easy access with radiation, 
because of its open geometry, and lends itself to form a two dimensional array on the 
same substrate. In addition, it has the advantage that well-known methods to produce and 
miniaturize it are available. 

As shown in Fig. Q it consists of a collection of circular electrodes printed on an isolating 
substrate. The simplest planar structure has two different electrodes to which voltages of 
opposite sign are applied. The trap is lying on the x — y plane and provides an electrostatic 
potential minimum along the z axis at a distance zq from the substrate. A strong homoge- 
neous static magnetic field of the order of few Tesla is aligned in the z direction and provides 
radial confinement. 




FIG. 2: Electrostatic potential with three electrodes {Vi = 2>V,V2 = — lOF, V3 = 2>V). One can 
see that the local potential minimum, in the ^-direction around p = 0, is almost quadrupolar even 
with the intrinsic asymmetry of the trap. 

An example of 3D potential is shown in Fig. |21for a configuration having three electrodes. 
The potential is obviously asymmetric and this implies that one needs to fine tune the applied 
voltages with more care than in a cylindrical trap. Another implication of its asymmetry is 
that one can move the position of the minimum along z. By varying the relative strength of 
the negative electrodes with respect to positive ones, a particle which is trapped around the 
minimum will vary its location under the control of the experimenter. This feature could 
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be useful to study the problem of decoherence caused by a metallic plate, which depends on 
the distance of the particle from the plate jl?! . 

The potential along the z axis can, in the case of adjacent electrodes, be obtained ana- 
litically from a Bessel-Fourier series expansion and reads 

=Y.y^\ ; - , , (1) 

where Ri-\ is the inner radius of the z-th electrode [Rn = 0, since the first electrode is a 

n 

disc). Vi is the voltage applied to electrode i and Nei^c. is the number of electrodes The 
minimum of that function cannot be evaluated analitically, but it can be shown numerically 
that zq^ the minimum position, has a value of the order of the radius of the central 
electrode. 

For electrons lying in the region near the trap axis, p -C the potential is well ap- 
proximated by Eq. ((H). Moreover, Eq. ((H) can be used to optimise the trap geometry and 
enhance its harmonicity. 



B. Harmonicity 

In an imperfect electrostatic field one is interested in cooling the electrons as much as 
possible so that they remain in the region near the minimum, where the potential is most 
resembling a quadrupolar one. A typical helium bath will thermalise electrons to a temper- 
ature of ~ 4 K and a dilution refrigerator at ~ 100 mK can drive them, on average, to the 
~ 200th energy level of their axial motion, for frequencies of ~ 10 MHz (which is the case 
when Vi ~ 0.01 V, -Ri ~ 1 mm). The latter temperature corresponds to an axial amplitude 
of ~ 20 /im and radial amplitude of the same order of magnitude. The amplitude of motion 
is thus much smaller than the characteristic trap size -Ri. Further cooling by pulses can be 
used in order to reduce the width of the axial oscillation. 

Detection of the electron axial motion can be performed by pickup of the induced image 
current in the central trap electrode, via a tuned resonance circuit of high quality factor 
Q. A Fourier transform analysis of the induced current shows a peak at the electron os- 
cillatory frequency, whose width AtUz is give n by the inverse time cooling constant due to 



the resistivity R of the detection electronics 
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With values for the tank circuit such as 



Q = 300, C = 7.5 pF and for ujz/2n = 10 MHz we have Aujz/2'n' ~ 1 kHz. Inside that peak, 



a dip is found, due to the electron which resonantly absorbs energy from the thermal noise 
in the circuit. Such a dip has a frequency width given by the axial anharmonicity. From 
numerical estimates we expect such a frequency width, for electrons which are thermalised 
to an environmental temperature of about 100 mK, to be, at least, two orders of magnitude 
narrower than that of the tank circuit |12i] . 

Moreover, to prevent detrimental effects on the computation due to electronic noise fed 
into the system, it is possible to decouple the axial motion of the electron from the detection 
circuit during gate operations. This is achieved by detuning either the external circuit or 
the axial oscillator. Only at the end of the computation, the axial frequency is brought 
into resonance with the detection circuit in order to perform the final read out of the qubit 
states. 

A superconducting solenoid provides the magnetic field along the z axis. A field strength 
of the order of few Tesla gives a cyclotron frequency of approximately 100 GHz which ensures 
that electrons will radiate their cyclotron energy, via synchrotron emission, in the time scale 
of 0.1 s (see Ref. Q]). Experimental observation of a trapped electron cooled down to its 
cyclotron ground state via radiation and in equilibrium with its cryogenic environment has 
been reported by Gabrielse and coworkers ^] under conditions similar to those mentioned 
here. 

Furthermore there is a way of coupling the spin motion to the axial one in order to know 
the spin state from a shift in the axial frequency. Such a coupling is realized when a quadratic 
magnetic gradient is applied. The method is well known and has been experimentally 
demonstrated for detection of the electron spin state in a hydrogen-like ion jlj]. In next 
sections we will show that a gradient, which is linear in the coordinates, can be used to 
provide an effective spin-spin coupling between different electrons. Therefore, if one desires 
to couple electrons and at the same time observe their spin state through measuring the 
axial frequency, both linear and quadratic gradients are needed. It is not a problem, in 
general, to tailor the required magnetic field. A ring made of magnetic material, such as Ni, 
has been already used for (yf-factor experiments. Such a ring provides a quadratic gradient, 
when seen from its symmetry centre. If the ring is displaced with respect to the trap axis, 
linear components will be seen by the trapped electrons in addition to quadratic ones. Other 
configurations are possible with the use of additional coils. 

With a typical nickel ring, the magnetic gradient is such that the electron axial resonance 
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frequency suffers a shift of approximately 10 Hz depending on the spin state. Considering 
the parameters for the tank circuit given above, a dip shift of 10 Hz inside a broad peak of 
1 kHz is easily detectable, furthermore when the absorption linewidth is 10 Hz - which is 
our case -. 




FIG. 3: Two dimensional array of planar traps consisting of printed electrodes on a base substrate 
made of isolating material. Each trap is supposed to be loaded with an electron cloud which is 
afterwards reduced to a single electron via a sequence of pulses. 



C. Array of planar traps 

The planar structure of the novel trap strongly suggests a 2D array of such traps, as 
schematically illustrated in Fig. El Thin-film technology can be used to place gold electrodes 
on an isolating surface with resolution much below the mm scale. A large number of traps 
can be embedded in a common isolating substrate and controlled electronically from the 
rear side. Ideally each trap is filled with a single electron, which interacts with neighboring 
particles via the Coulomb force. This provides a natural multiparticle scenario similar to 
the case of ions in a linear Paul trap, with the advantage of controlling each trap parameter 
independently (inter-particle distance, coordination number, electron resonance frequencies, 
. . . ). With the planar trap this is implemented in a very straightforward way. 

The substrate shown in Fig. 01 has the size of a coin and each trap a dimension of order 
1 mm. This size gives an axial frequency up to 500 MHz for applied voltages of a few 
Volts. The space between traps has to be filled with a common grounded electrode, so as 
to isolate the potentials of each trap from the others and also to prevent charging up of 
the isolating substrate. Thin-layer techniques allows to produce electrodes which are almost 
monocristalline, and so to reduce significantly the decoherence effect due to patch fields. 
On the other hand, quantitative measurements of decoherence caused by proximity to the 
electrodes could be performed by moving the electron along the z axis (see jl2i]). 
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Currently, we are testing a prototype of planar trap. The idea is to have only one trap, 
with a radius of 2 cm, printed with thick-layer technology over a AI2O3 substrate. A magnetic 
field of 100 G is foreseen. This preliminary trap will try to demonstrate confinement and 
also the possibility to excite each degree of freedom of the electrons. Figure E] shows the 
design of such a test trap. It can be seen that there is a split electrode, which will provide 
quadrupolar excitation. In a further step, we plan to build a miniaturised array of several 



FIG. 4: Test trap with total diameter of 35 mm. Black regions denote electrodes and white regions 
denote isolating surface. As said before the isolating surface has been mostly covered by electrodes 
to avoid charging up. This configuration has three active electrodes plus an external grounded one. 
The red circle in the center is a hole through which electrons are loaded from the rear side. 

traps (three for example) in a cryogenic environment of about 100 mK with a magnetic field 
of 7 Tesla. The traps will have a diameter of around 0.5 mm. The biggest challenges will 
be those of coherent and accurate control over each electron plus a sufficient supression of 
all sources of decoherence, both haunting every experiment in quantum computation. 

III. BUILDING AN ARTIFICIAL MOLECULE 
A. The magnetic gradient 

Let us consider a planar trap at a distance Xq from the center of the substrate along the x- 
direction. We choose the origin of our reference frame at the substrate center and the 2;-axis 
orthogonal to the substrate plane (see Fig. E)). Now we suppose to add the inhomogeneous 
magnetic field 





(2) 
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FIG. 5: Schematic draw of a planar trap at the distance xq from the center of the substrate. 



This field produces a linear magnetic gradient and has rotational symmetry with respect to 
the z-axis. 

The total magnetic field acting on the trapped particle is B = Bq + Bi where Bq is the 
uniform confining field directed along the 2;-axis. Hence, the total field B acting on the trap 
center is the sum of the uniform field Bq along the z-direction and the magnetic gradient 
—6x0/2 along the x-direction. Consequently B, at the trap center, has modulus Be = 



Bq + (62xq)/4 and direction forming an angle with the z-axis. Note that Be depends on 
the distance |xo| of the trap from the substrate center. Now we rotate the reference frame 
around the y-axis so that the new z-axis corresponds to the direction of the total magnetic 
field B at the trap center. With this rotation, if 6|xo|/2i?o ^ 1, the total field with respect 
to the new coordinates can be written, in good approximation, as 



B 



[Be + bz)k 1 - yj 



(3) 



The vector potential of the total magnetic field B applied to the electron is 



A ~ ^ (5c + bz) [{x - xo)j - yi] . (4) 

We suppose to work with a uniform field Bq of few Tesla, a magnetic gradient b of about 
50 T/m and a trap substrate with length of the order of 10~^m. Consequently we have 
b\xo\/2Bo < 10-\ 

Let us write the Hamiltonian of the trapped electron. By taking into account the trapping 
potential we have 

H= ^^-'^^\ eV-^. .B, (5) 
11 



where V, rrie, e, g, and cXi are, respectively, the trapping potential, electron mass, electron 
charge, giromagnetic factor, and Pauli matrices. In the limit b\xo\/2BQ ^ 1 we can neglect 
the changes in the quadrupole potential form due to the rotation of the reference frame 
around the z-axis and write 

2 {x-xoY + y'^ 



T/ ^0 
y ~ — 



(6) 



We can define the axial frequency ujz = y2eVQ/ (rrieP), in terms of the applied potential 
difference Vq and of the characteristic trap size i, and recast the z-part of the spatial Hamil- 
tonian of the electron as 



The presence in the Penning trap of the magnetic gradient along the z-direction makes 
the cyclotron frequency depend on the particle 2;-position. Indeed, we define the cyclotron 
frequency as 



\e\{B, + bz) 



me 

and introduce respectively the cyclotron and magnetron ladder operators 

1 I ImCjc . / 2 . 

^[(x-xo)-zy] + W^(p, + zp.) 



1 / muOr 



[{x - xo) + iy] 



— {Py - iPx ) 



(8) 

(9) 
(10) 



obeying the commutation relation [oj, a]] = with j = c, m. The frequency uOc is defined 



as 



- 2ul 



(11) 



so it depends on the z-coordinate too. 

The part of the spatial Hamiltonian of the electron involving x-y coordinates is 



H ~ 



1 



2mp 



e{B^ + bz) 
Px + y 



1 2 



+ 



2mp 



e{Bc + hz) 
Py [x-xo) 



n 2 



-meujl[{x - xof + y\ 



(12) 



By using the ladder operators, Eqs. Q and (fTUj) . we can write 



Hxy = -huJm ( al^am + 2 j + (^«I«c + 



(13) 
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where we have introduced, respectively, the magnetron and cyclotron frequencies 

ujm{z) = , (14) 

= 2 ' ^ 

both depending, by means of relations (jH)) and on the z-coordinate. Notice that we 
indicate the explicit dependence of the frequencies Uc, ojc, ^'ci ^i^d Um on the coordinate z 
only in the definition formula. 

The uniform magnetic field of few Tesla, we suppose to apply, gives a cyclotron frequency 
at the trap center of the order of 100 GHz. This choice permits to have, after thermalization 
with a trap environment at about 100 mK, the cyclotron motion in the lowest energy state 
[2^. Furthermore we choose the axial frequency in the range of MHz, giving a magnetron 
frequency at the trap center of the order of hundreds of Hz. 

The spatial Hamiltonian, Eq. (fT^ . is formally equivalent to the corresponding part of 
the Hamiltonian describing an electron in a conventional Penning trap, without any mag- 
netic gradient 0]. However, we point out that, with respect to the case of a trap with- 
out magnetic gradient, there is a dependence of the magnetron and cyclotron frequencies 
on the 2;-coordinate. Hence, we have a coupling between the axial motion and the cy- 
clotron/magnetron motions. Furthermore, as we shall see, the magnetic gradient introduces 
an interaction between the spatial motion and the spin motion of the electron. This coupling 
between the external and internal degrees of freedom becomes evident by considering the 
part of the electron Hamiltonian involving the spin motion 

geU gh gfi\e\b 

Hs ^ cr ■ B = —uJc(T^ - — {a^x + cjyy). (16) 

4me 4 8me 

We generally suppose, in our system, the axialization of the electron motion. This condition 



22 1 . The motion of 



has been experimentally obtained with ions confined in a Penning trap 
the electron is axialized when its amplitude in the x-y direction is much smaller than that 
in the 2;-direction. The axialization permits to neglect, in the above Hamiltonian, the term 
proportional to + ayy. Indeed, with the typical values chosen for h and Bq, the spin 
state transition probability due to this term is negligible in the axialization condition. This 
can be proved by applying time- dependent perturbation theory. 

Hence, the global Hamiltonian of the electron can be rewritten as 

13 



where the frequencies Um, ^'^ and lUc, as defined in Eqs. (|HI), (HH), (d), and (dS)), depend 
on the z-coordinate. 

We suppose that, throughout the electron motion, the conditions <^ cUc and b\z\/2BQ <c: 
1 are verified. Hence, by posing uj'^ ~ uJc, — ^l/'^^c and expanding the cyclotron and 
magnetron frequencies in terms of the coordinate z, we write the above Hamiltonian as 

Ti 

+ huj^e{a, + 4) ( -^aj„a„ + aja^ + f ) , (18) 
where we have introduced the ladder operator 



rripUJ 



^ + hU^-—P^. (19) 



(20) 



2h V 2hmeUi 

the parameter 

^ \e\b I h 

nieUJz V '^rUeU^z ' 

and the frequencies Uco = (|e|-Bc)/'"^e, i^mo = ^^'i/'^^co, ^^so = (5"^co)/2. 

In Hamiltonian (fTH|) we have also redefined the equilibrium position of the electron in the 
2;-direction, as a consequence of the shifts due to the zero point energy of the cyclotron and 
magnetron oscillators. 

The parameter e gives, substantially, the ratio between the change in the electron spin 
energy due to the magnetic gradient and the trapping energy in the z-direction. Indeed we 
can write 

hd;,uJsAz 



where 

2m 



dzUJs = t; (22) 



and 



Az=J—^ (23) 

V 2meUJz 

is the amplitude of the electron axial motion in the ground state. Hence dzUJsAz is roughly 
the spin frequency variation due to the magnetic gradient when the electron axial motion is 
in the ground state. 

Note that the frequencies lJcq, iVmo and Uso depend on the position of the planar trap on 
the substrate. Indeed they depend on the distance xq of the trap center from the symmetry 
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axis of the magnetic gradient. For instance, when the magnetic gradient 6 ^ 50 T/m, 
two electrons in neighboring traps, separated by a distance of the order of 10^'^ m, are 
characterized by spin resonance frequencies that differ from each other of few MHz. This 
value, though small in comparison with the typical spin frequency uJsq/2'k ~ 100 GHz of a 
single electron, is enough to individually address the spin qubits via microwave radiation. 
Moreover, the same substrate can accomodate tens of qubits, with frequencies spread over 
a range of hundred MHz. 

If the electron motion is axialized and the cyclotron oscillator remains always in 
the ground state, we can neglect, in Hamiltonian ()18|1 . the term proportional to (a^ + 
a\)[{umo/^cQ)<An'^m + <i\cLc\. Hcucc, Under these conditions, we can write the Hamiltonian of 
the system as 

H ~ -/iWmoaL^m + ^LUcoaJac + ^w^ajaz + -c^soc^z + ^^^^2^(0^ + 4)'^^? (24) 
showing the coupling between the axial and the spin degrees of freedom. 



B. Effective spin-spin coupling 



We now move to a linear array of electrons confined in planar Penning traps along the 
X-axis. If we add a linear magnetic gradient, the Hamiltonian of the system can be written 

as 



H 



N 

E 

1=1 



[pi - eAi) 



geh 



N 



+y — - 



(25) 



2me 4me 
where the subscript i refers to the i-th electron of the array. 

Hence the magnetic field Bi, with vector potential Ai, acts on the electron trapped in the 
i-th site. This field is the sum of the uniform confining field orthogonal to the planar trap 
substrate and the inhomogeneous field producing a linear magnetic gradient [see derivation 
of Eq. Q in the previous subsection] 



Bi 



(26) 



In the above equation -Bc.ik is the total field Bi at the center of the i-th trap and Xi^ is 
the x-coordinate of the center of the i-th trap. The origin of our reference frame is at the 
center of the trap substrate and the 2;-axis is the symmetry axis of the magnetic gradient. 
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We have i?c,i = y-Bo + {b'^^ifl)/'^ where \xifi\ is the distance between the symmetry axis of 
the magnetic gradient and the center of the i-th site. In writing Eq. ()26|) we assumed that 
b\xo,i\/2Bo < 1 for any i = 1, 2, ... , A^. 

Similarly to the single electron case, we also have 

Ai ~ + bzi)[{xi - Xi,o)j - yd]- (27) 



and 



^2 _ i^i ~ ^ifi) + Hi 



(28) 

Now, as already done for a single electron in a planar trap in the presence of a magnetic 



gradient, we define the frequencies ujc,i{zi) = |e|(-Bc,i + bzi)/me and ujc,iizi) = yi^c,i ~^ ^cu^ 
and introduce the operators 



a™. = 2 j y ^ [(^. - ^..o) + WA - ^ ^— (p..- - ^P..) J , (30) 



a. 



rripUJz . / 1 



obeying the commutation relation [aj,aj] = with z,j = c,m,z. 

If we indicate with if^*^ the part of Hamiltonian ()25|) not including the Coulomb inter- 
action, we can write 

N / n ^ 



^ 9 



+ 7^^^^(«2,i + (32) 
i=i ^ 

where we have introduced the frequencies cjco.i = (|e|-Bc,i)/'^e, ^^mo,i = ^1/'2^co,i, ^so,i ^ 
{gUco,i)/2 and assumed that, throughout the electron motion, the conditions <^ ujc^i and 
b\zi\/2Bo <^ 1 are verified for any i = 1,2, . . . , N. The parameter e has been defined in Eq. 
()20|) . In writing Hamiltonian ()32|) we have also supposed that all the electron motions are 
axialized and that the cyclotron oscillator remains always in the ground state. 

Let us consider the part of Hamiltonian ()25|1 involving the Coulomb interaction between 
the two electrons at the sites i and j. By indicating it with Hfj we have 

= , (33) 

Aneo^ixi - xj)^ + (jji - yj^ + {zi - zj^ 
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which we can recast as 



1 + ^ ^ + ^ ^ + 



2 (z, -zA^ 



d - ■ S ■ ■ d^ ■ 

where Axi = Xi — Xi^o and (ijj = \xifi — Xjfi\. 

If the oscillation amplitude of the two electrons is much smaller than the average sepa- 
ration dij between them, we can expand the interaction Hamiltonian, Eq. in a power 
series and retain the terms up to the second order 

2 2 

Furthermore, if we suppose that the electron motions are axialized, we have 

2 2 

47reoafj oneodlj 

The first term in Hamiltonian (|36p gives a displacement of the equilibrium position of the 
electrons along the x-axis. This small shift is of order of e'^ / (ineod'^jmeUj'^Q i) and is a 
consequence of the Coulomb repulsion. Hence we can effectively remove the first term 
in Hamiltonian ()36|) by redefining the centers of the two traps. We also recall that, actually, 
this effect is more pronounced for electrons placed at the extremities of the array. Indeed, 
as a consequence of the symmetry of the system, particles trapped near the center of the 
array undergo much smaller shifts. 

The second term in Hamiltonian (j36p . involving the z-coordinate of the electrons, repre- 
sents an effective dipole-dipole interaction between the i-th and j-th electrons. By devel- 
oping the square we have terms proportional to zf and z|. They produce small shifts on 
the axial frequencies of the two electrons. Therefore we can take into account these terms 
by appropriately redefining the axial frequencies of the two electrons. The remaining term, 
proportional to ZiZj represents the Coulomb coupling between the axial motion of the two 
electrons. After having appropriately redefined the trap centers and frequencies we can write 

g2 

H^j - 4^g^^3 = ^^i,Mz,i + al,i)i.(^z,j + 4j), (37) 
where ^ 

~ 47reo(iij (ciij) ' ^^^^ 
with Az being the ground state amphtude of the axial oscillator, introduced in Eq, 
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Now we perform, on the global Hamiltonian of the system H = H^^ + I]i> j Hf- , the 
unitary transformation H' = e^He^^ [5] with 

^ 9 

S = Y1 - (^z,i)(^z,i- (39) 

The transformed axial operators are 

az,i a^^i - e^(Jz,i- (40) 

Let us consider the transformed part of the Hamiltonian not including the Coulomb 
terms. It can be written, after dropping constant terms, as 



(41) 



Note that, with the unitary transformation, we have formally removed the interaction be- 
tween the axial and spin motions. Differently the Coulomb terms in Hamiltonian ()25|) 
involving the electrons i and j transform as 

From the above Hamiltonian we see that the unitary transformation produces terms of the 
form hlg"^ /A)e'^C,i,jCrz,iO'z,j, which represent an effective coupling between the spin motion of 
different particles. 

By applying another unitary transformation, similar to Eq. (I39j) . also the extra terms 
in Hamiltonian (j^ . proportional to az^iia^j + ^Ij), result in additional spin-spin coupling 
terms. As a consequence we would have corrections to the coupling strength of the spin-spin 
interaction the order of he^S^f j/ojl- However, we assume ^ cu^, so that the effect due to 
the terms proportional to (yz,i{ciz,j + ,,•) is negligible. 

Since in our scheme quantum information is encoded only in the spin motion of the 
particles, we do not consider the spatial part of the transformed Hamiltonian, but just the 
spin-part, given by 

K - Y.2^so,i(Tz,i + Y.2'^Ji,i^z,t(Jz,j- (43) 

i=l i>j 

Note that the above Hamiltonian is analogous to the nuclear spin Hamiltonian of the 
molecules used to perform NMR quantum computation Consequently, with our system, 

18 



we can implement universal quantum processing by using the same techniques developed in 
NMR experiments (see next section). 

In Hamiltonian (|43p the spin frequencies are 



and the coupling constants, defined as in NMR experiments, are 



A. = 7^eM^ =7^7^7T^h^ ^^T^- (45) 



Hence, by changing specific system parameters, fully under the control of the experimenter, 
we can adjust both the spin frequencies and the coupling constants. Indeed, according to 
the above relations, the spin frequencies of the particles and, consequently, their detunings 
depend on the uniform magnetic field intensity, the magnetic gradient strength and the 
distance of the electrons from the substrate center. Similarly, the coupling constants can be 
modified by changing the gradient strength, the axial frequency and the distance between the 
particles. We also note that the coupling constants are proportional to l/rf^^. Hence, if the 
particles are equally spaced in the array, we obtain a uniform coupling strength for nearest- 
neighbor electrons, while the spin-spin interaction decreases rapidly with the increase of their 
distances. For instance, the coupling strength between two electrons, that are a distance Id 
apart from each other, is just 1/8 of the one for nearest-neighbor electrons. In the considered 
dipole limit, we achieve the same uniform interaction strength between neighboring electrons 
and reduction for non-nearest-neighbor J couplings of Ref. j^, where, however, this result 
requires a careful adjustment of both the inter- ion separation and end-trap strength. 



IV. UNIVERSAL QUANTUM LOGIC GATES 

In this section we describe how to implement, in our system, universal quantum compu- 
tation. 

As mentioned in the previous section, Hamiltonian is substantially analogous to the 
one describing NMR molecules |lj . The only difference is the fact that in the NMR systems 
we have nuclear spins instead of electron spins. Therefore, universal quantum processing 
with trapped electrons, in the presence of a magnetic gradient, can be performed by using 
techniques similar to those already developed in NMR experiments 
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In particular, Hamiltonian ()43|) represents a set of electron spins, each one having a 
different precession frequency, which interact by mutual couplings. If we encode a qubit in 
the spin degree of freedom of each particle, we can use electromagnetic pulses with appro- 
priate frequencies to selectively manipulate the information stored in the spin state of each 
electron. Hence, the detuning between the different electron frequencies plays the same role 
of the chemical shift in NMR molecules. 

Universal two-qubit gates, in our system, are achieved by means of the mutual coupling 
between the electron spins. Indeed, this interaction has the same form of the J-coupling 
between nuclear spins in NMR molecules and can be used in a similar way to perform 
two-qubit gates. 

The set of unitary transformations consisting of single-qubit gates plus C-NOT gates is 
computationally universal. Let us describe, in detail, how to perform, in our system, these 
operations. If we apply a small transverse oscillating magnetic field resonant with the spin 
precession frequency Usoj of the j-th electron 

Bd(t) = Bd[icos{ujso,jt + 6) + J sin(cj,o,it + 0)], (46) 

the relevant part of the system Hamiltonian becomes, in interaction picture with respect to 
the unperturbed Hamiltonian Hso = Z]ili(^/2)co'so,i'^z,i ^^nd in rotating wave approximation, 

Hir^ ^ (a+,e-^^ + a_,e^^) , (47) 

where x = g\e\Bd/{2me) and a±j = {a^j ± icryj)/2. 

In the derivation of the above Hamiltonian we neglected the spin-spin coupling terms present 
in Hamiltonian ()43|) . since we assume Jij -C x for hi = 1, 2, . . . , A^. Furthermore we 
suppose that the oscillating field is so small that ^ foc/j j for any i, j = 1, 2, . . . , A^. This 
last condition, giving a Rabi frequency x much smaller than the detuning between the spin 
frequencies of the z-th and j-th electrons, allows for the selective frequency addressing of 
each particle. 

If the small transverse magnetic field is applied for a time t, it produces a rotation on 
the spin state of the j-th particle 

i).-^e-^m(f)|T)„ (48) 
T),-ze^^in(^)|i),. (49) 
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i). - cos (I) 
T). - COS (I) 



It can be shown that with an appropriate combination of these operations, one can perform 
any single-qubit gate on the spin qubit of the j-th electron. We define the interaction 
produced by the Hamihonian (jTTj) . apphed for a time t, as a Psjix^i ^) pulse. Hence, we can 
perform single-qubit gates on each electron, in a selective way, by appropriately tuning the 
frequency of the oscillating magnetic field. 

Let us now consider the implementation of C-NOT gates. A natural way to perform 
this two-qubit gate is to use a three gate circuit. In particular, we apply in sequence, an 
Hadamard gate on the target qubit, a controlled vr phase-shift gate followed by another 
Hadamard gate on the target qubit. However, it is preferable to avoid Hadamard gates, as 
they are difficult to implement. Hence, the Hadamard gates are conveniently replaced by an 
inverse pseudo-Hadamard gate and a pseudo-Hadamard gate . These two gates are realized 
by applying respectively a ps(7r/2, 7r/2) pulse and a ps{7f/2, — vr/2) pulse. The controlled vr 
phase-shift gate can be achieved, in different ways, by applying appropriate sequences of Ps 
pulses. For example, in a two-qubit system, one possible sequence for implementing this 
gate consists of four ps pulses and two appropriately timed periods of free evolution In 
particular, we apply: a free evolution of the system for a time of l/4Ji^2, a Ps{tt, 0) pulse, 
another free evolution of the system for a time of 1/4 Ji 2, a Ps{7t/2, 0) pulse, a Ps(7r/2, 7r/2) 
pulse and a ps(7r/2,0) pulse. All the above pulses should be applied to both electron spins 
and the application time of each p^ pulse should be much smaller than the time duration of 
each free evolution period. 

With systems having more than two qubits, the efficient implementation of C-NOT gates 
requires the application of more complicated sequences of operations. This is necessary in 
order to avoid errors due to the couplings with electrons not involved in the gates. Indeed, 
from Hamiltonian ()43|) we see that there are mutual couplings among all the particles of the 
system. However, the sequence required to efficiently implement the C-NOT gates, when 
we have more than two qubits, consists only of Po pu lses acting on specific particles and 
appropriately timed periods of free evolution j22l, |2J|. In NMR context this technique is 
known as refocusing. 

We also remark that in our system the coupling strength decreases rapidly with the 
increase of the inter-particle distance. Indeed from equation we have Jij oc ^/d^j- 
This fact permits to simplify the refocusing sequences since the interaction between distant 
electrons can be neglected recall that in this case the logic gates between distant 
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electrons can be performed by means of swap gates, which move the quantum information 
among the particles. A swap gate is reahzed by applying three C-NOT gates, where the two 
qubits play alternatively the role of target and controller. 

Let us give realistic estimates for the values of the electron detunings and spin-spin 
couplings achievable in our system. We consider a linear array of 10 electrons with inter- 
particle distance of the order of 1mm. We suppose to apply a uniform magnetic field of few 
Tesla, giving ojsjl'K ^ 100 GHz, a magnetic gradient of about 50 T/m and assume an axial 
frequency uJz/2t: ~ 10 MHz. In these conditions we obtain a frequency detuning between 
neighboring particles of few MHz and a spin-spin coupling with strength J fa 20 Hz. These 
values for the detuning and the coupling strength are of the same order of the corresponding 
quantities in NMR systems. 

V. CONCLUSION 

A novel concept of open planar Penning trap makes it possible to design and build-up 
a scalable system for quantum computation, consisting of trapped electrons in vacuum. 
Single particles are confined in a linear array of traps, deposited on the same substrate. 
A magnetic field gradient across the traps allows for frequency addressing of each qubit, 
stored into the electron spin as in NMR quantum computers. Moreover, the magnetic field 
gradient couples internal (spin) and external (motional) degrees of freedom of the same 
particle. Thanks to the Coulomb interaction among the charged particles, this coupling 
effectively amounts to a direct spin-spin interaction, with tunable coupling strength. In 
the limit of relatively large inter-electron spacing and small axial oscillation amplitude, we 
obtain an analytical expression for the J coupling, which allows for an immediate evaluation 
of the interaction strength in terms of the relevant system parameters (trap distance, axial 
frequency, and applied magnetic gradient) . We emphasize that the J coupling is proportional 
to where d is the inter-trap distance, thus greatly reducing the interaction between 

non-nearest-neighbor electrons. 

This way, the system of singly trapped electrons is formally identical to a NMR molecule, 
suitable for quantum information processing. Hence, the well established and developed 
techniques of NMR spectroscopy can be extended and applied to our system. Qubit ma- 
nipulation is achieved via appropriate sequences of microwave pulses. The qubit read-out 
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is performed either by axial frequency detection, as in conventional Pennning traps, or by 
capacity and charge measurements as in semiconductor quantum dots. 

The advantages over NMR systems are obvious: the number of qubits is not limited 
neither by the molecule size - the same substrate can easily accomodate a large number 
of microtraps - nor by the frequency range - the typical spin resonance frequency lies 
in the GHz, whereas the detuning between neighboring particles is of the order of a few 
kHz -, and the system is truly quantum and, at the same time, well isolated from the 
environment. In addition, the geometry of the system can be designed and optmized at will, 
using more complicated two dimensional arrays, with possible applications to the study 
and the simulation of quantum systems like the Ising model. Finally, Penning traps could 
be loaded with protons as well, allowing the same NMR spectroscopy experiments with 
hydrogen, but in a much more controlled and clean environment. 
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